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EXPERIMENTAL FLUTTER STUBI OF A 


WING-FUSELAGE CONFIGURATION AT A MACH NUMBER OF 15.4 

AND COMPARISON WITH THEORY 

By Robert C. Goetz and John L. Sewall 
Langley Research Center 


SUMMARY 



An investigation was conducted in helium flow at a Mach number of 15.4 to 
provide hypersonic flutter data on a low-aspect-ratio, flexible, wing-fuselage 
configuration. The models were sting-mounted in the 24-inch-diameter leg of 
the Langley Mach 15 hypersonic aeroelasticity tunnel on a support system which 
simulated the free-flight condition by permitting symmetric freedoms of vertical 
translation and pitch. The fuselage was cylindrical with a semirigid spheri- 
cally blunted nose-cone section, and the wing was a 73° clipped delta with a 
modified trailing edge and a slab airfoil section. The model geometry and mass 
distribution were held constant while the ratio of wing stiffness to fuselage 
stiffness was varied. 

The results indicate that increasing the wing-to-fuselage frequency ratio 
from 0 to about 1.0 is destabilizing, whereas a further increase in frequency 
ratio is stabilizing. However, the flexible model never becomes as stable as 
the rigid-fuselage model over the test range of wing-to-fuselage frequency 
ratios. From these results the fuselage flexibility was concluded to be an 
important flutter parameter. The stiffness of the model support system was 
at a level low enough to have no measurable effects on the flutter character- 
istics of the models in the present investigation. 

Flutter calculations for some of the models were performed by application 
of a flutter analysis based on the first four coupled symmetric modes and on 
quasi-steady aerodynamic theory for approximating aerodynamic forces acting on 
wing and fuselage. Calculated flutter-speed trends were similar to experi- 
mental trends but were higher by as much as about 80 percent . The lack of a 
sufficient number of flexible wing-fuselage modes prevents a true assessment of 
the aerodynamic approximations used in the analysis. The use of a combination 
of Newtonian theory for the fuselage nose cone and piston theory for the rest 
of the configuration resulted in essentially the same flutter solution as that 
obtained by using piston theory for the entire wing-fuselage combination, a 
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INTRODUCTION 


Methods for treating the flutter characteristics of wing-fuselage combina- 
tions at high speeds have been formulated analytically in references 1 to 3j 
however, corroborating experimental flutter studies of similar wing-fuselage 
combinations are very scarce. Most experimental hypersonic flutter studies 
(refs. 4 to 8) have been devoted to specific wing-shape effects such as airfoil 
shape, leading-edge sweep, thickness- chord ratio, and leading-edge bluntness. 
These studies have a limited application to designs of proposed hypersonic 
vehicles in which the fuselage and the wing are merged so that it is impractical 
to treat them as separate components. In such designs the complete wing- 
fuselage combination must be considered. 

Accordingly, the present investigation was undertaken to provide experi- 
mental hypersonic-flutter data on a low-aspect-ratio, flexible, wing-fuselage 
configuration. The purpose of this investigation was to assess the importance 
of fuselage flexibility on the flutter characteristics of the wing-fuselage 
combination. The model configuration, tested in helium flow at a Mach number 
of 15.4, consisted of a body of revolution, or fuselage, with a centrally 
mounted wing. The fuselage was cylindrical with a spherically blunted nose- 
cone section, and the wing was a 73° clipped delta with a modified trailing 
edge and a slab airfoil section. Models of this configuration were sting- 
mounted on a support system that simulated the free-flight condition by per- 
mitting symmetric freedoms of vertical translation and pitch. The primary 
variable of this investigation was the wing-to- fuselage stiffness ratio. 

A coupled-mode flutter analysis based on quasi- steady aerodynamic theory 
was applied to models for which measured mode shapes were available, and com- 
parisons are made between calculated and measured flutter results. Both piston 
theory and Newtonian theory were used to approximate the oscillating aerody- 
namic forces acting on wing and fuselage, and wing-fuselage aerodynamic inter- 
ference effects were neglected. 


SYMBOLS 


A ij 
a 

bo 
C A 

Cl = oCl, 


aerodynamic coefficient in flutter analysis (see eqs. (A3) and (A 5 )) 
speed of sound, feet per second 

wing reference semichord, 0.642 foot for all models (see fig. l) 

Axial force 


axial-force coefficient, 


<3*M £ 


•a 


c l, 


'a 


lift- curve slope 


2 


c l 

C m 

C N 

F 

f 

g 

h i,hj 

3 

k 

Z M 

Z F 

^X 

M 

Mi 

mt 

mp 

my 


rolling-moment coefficient, ^ ^ ^ mQme 

<% r c 

pitching-moment coefficient about nose of model, momen t 

normal-force coefficient, formal force 

dZ M 

series defined by equation (A4) 
frequency, cycles per second 

2jt 

structural damping coefficient in flutter eigenvalue (see section of 
appendix following eq. (Al)) 

normalized mode shape for ith and jth mode, respectively 
integers identifying modes used in analysis 
reduced frequency, 

reference length of model, 1-75 feet for all models 
fuselage length, 2 x 7p 

reference length in x-direction, 2b 0 in the present paper 

length of fuselage nose cone including rounded nose, Z x 7^, 

0.521 foot for all models 

Mach number 

generalized mass in ith mode (see eq. (A2)) 

total mass, slugs 

mass of fuselage per unit length 

mass of wing per unit area 

dynamic pressure 

radius of fuselage at any chordwise station, l x r 


3 



radius of cylindrical fuselage section, 7 x r c , 0*125 foot for all 
models 

wing semi span measured from fuselage center line 
local airfoil thickness 
velocity, feet per second 

chordwise coordinate, positive rearward, measured from fuselage nose, 

V 

spanwise coordinate, measured from model center line, sy 
angle of attack (positive nose up), degrees 

ratio of specific heats, 5/5 for helium (see eqs. (A3) and (A4)) 
semivertex angle of fuselage nose cone, 12 degrees 

angle between local tangent to fuselage surface and longitudinal axis 
of fuselage 

) 

— (see section of appendix following eq. (Al)) 

nn 

mass-ratio parameter, — 

P^IT 

density, slugs per cubic foot 

complex flutter eigenvalue, (l + ig)> (see section of appendix 

following eq. (Al)) 

angular frequency, 2rtf, radians per second 
reference frequency, radians per second 

flexible-to-rigid wing-mounting stiffness ratio, or wing-flapping 
stiffness ratio 

wing-to-fuselage frequency ratio 


wing-to-fuselage stiffness ratio 


V 


flutter- speed- index parameter 


Subscripts: 

ac aerodynamic center, measured from model nose, fraction of model 

length 

eg center of gravity, measured from nose, fraction of model length 

f pertaining to flutter mode 

s pertaining to wing mounted on rigid support 

W pertaining to uncoupled wing flapping mode 

A pertaining to wings mounted on body with bending stiffness A 

B pertaining to wings mounted on body with bending stiffness B 

C pertaining to wings mounted on body with bending stiffness C 

A bar over a symbol denotes a nondimensional quantity. 

EXPERIMENTAL PROGRAM 


Description of Models 

The model configuration tested in the present investigation consisted of 
a body of revolution with a centrally mounted wing which had a rounded leading 
edge swept 73 ° and a slab airfoil section. The model geometry is shown in fig- 
ure 1. The fuselage consisted of a semirigid spherically blunted nose- cone 
section which had integral longitudinal inserts along the wing plane. These 
inserts extended back from the nose-cone section to constitute the main longi- 
tudinal bending stiffness of the flexible portion of the fuselage. A coil 
spring was wound through, and brazed to, the inserts in order to form a cylin- 
drical body without contributing significant additional longitudinal bending 
stiffness. The coil spring was covered with a layer of rubber (dental dam) in 
order to provide a relatively smooth aerodynamic surface. In figures 2(a) to 
2(d) the model can be seen in the various stages of construction. 

The wing planform was a 73 ° clipped delta with a modified trailing edge. 

The wing was constructed of balsa wood with two steel flexure-beam inserts 
oriented in the spanwise direction. The exposed portion of flexure beams was 
3/l6 inch wide and 1/8 inch long, and the thickness was varied in order to 
change the wing flapping stiffness. These flexure beams were bolted to the 
integral longitudinal inserts at two fuselage stations. (See fig. 2(e).) Model 
geometry and mass distribution were held constant while the ratio of wing stiff- 
ness to fuselage stiffness was varied. A rigid-fuselage model was constructed. 
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with the wings connected in the same manner as they were on the flexible model, 
in order to assess the effect of fuselage flexibility on the flutter character- 
istics. Figure 2(f) shows the rigid- fuselage model with the upper portion of 
the cylindrical fuselage rotated out of place. 

Each model is designated by a capital letter, an integer, and a lower case 
letter. The capital letter indicates the level of the fuselage longitudinal 
bending stiffness: A is stiff er than B, and C is rigid. The integer indi- 
cates the level of wing flapping stiffness: 5>k>3>2>l. The lower-case 

letter indicates the level of model support stiffness: a > b > c. For example, 
C-l-c means that the fuselage is rigid, the wing flapping stiffness is the 
weakest, and the model support system is the weakest. 


Model Support System 

The models were sting mounted in the tunnel on a flexure-pivot support 
system. The support system consisted of a pair of parallel flexure springs 
clamped on the sting, with a cross-pivot spring connected to the free ends. The 
free end of the cross-pivot spring was connected to the model by a vertical 
rigid rectangular shaft. Figure 3 shows a view of the model support system 
partially assembled. Three levels of support -system stiffness were used in 
order to determine whether the support system was in any way influencing, or 
contributing to, the model flutter characteristics. At all three levels the 
system had a translation-to-pitch frequency ratio of about 2/3* 


Physical Properties 

The total mass (excluding the mass of the model support system) and natural 
frequencies of the various models are listed in table 1(a). The weight distri- 
bution of the model support system is shown in figure 4(a). Figure 4(b) shows 
the weight distribution, obtained by weighing the various model components, for 
a typical flexible fuselage. A typical measured wing-weight distribution is 
presented in figure 4(c); this distribution was determined by cutting the wing 
chordwise and spanwise as indicated in the numbered pattern in figure 5 and 
weighing each segment. It should be noted that figure 4(c) also includes a 
table of the concentrated mass for each wing due to the steel flexure beams 
which extend over the entire wing span at their appropriate chordwise stations. 
The difference between the total mass of the rigid- and flexible-body models as 
shown in table 1(a) was in the cylindrical portion of the fuselage, and was 
uniformly distributed. This difference changes the location of the model cen- 
ter of gravity as is shown in figure 5* The cross-pivot location (pitch axis) 
of the support system was 0.4152^ (measured from the nose) for all models. 

All models were vibrated on the model support system with an interrupted- 
air- jet shaker to determine their natural frequencies and mode shapes. Mode 
shapes were measured by using a pivoted- coil-galvanometer element whose output 
was recorded on an oscillograph while the models were vibrating with constant 
amplitude in their natural modes. These displacement measurements were taken 
at various stations on the model (circled coordinates in fig. 5); and the 
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results were normalized with respect to maximum displacement in a given mode. 
These results for several of the models are presented in figures 6 to 9* A 
brief summary of the model motion in each given mode follows: 

(1) First mode (model translation) - The model behaved essentially as a 
rigid body with vertical translation predominating. 

(2) Second mode (model pitch) - The model behaved essentially as a rigid 
body with pitching predominating. 

(3) Third mode (model rolling) - The model oscillated about the model 
center line in an antisymmetric fashion. This antisymmetric mode has not been 
included in figures 6 to 9 because it was not regarded as contributing to the 
symmetric flutter motion. 

(4) Fourth mode (fuselage bending) - The fuselage had a definite deforma- 
tion along its length, and the wing motion was symmetric but out of phase with 
the fuselage motion. (The omission of this mode in fig. 8 is due to the rigid- 
ity of the fuselage for model C-l-a.) 

(5) Fifth mode (wing-fuselage combination) - The maximum displacement 
occurs at the wing tip. The wing motion is symmetric but out of phase with the 
fuselage motion. The wing is relatively rigid over its span, whereas the fuse- 
lage is flexible. 

In addition to the foregoing coupled modes, an uncoupled mode was meas- 
ured in order to investigate the effect of wing flapping stiffness. The fuse- 
lage was restrained along two planes tangent to its upper and lower cylindrical 
surfaces, while the wing flapping displacement was measured. The frequency of 
this mode was lower than that obtained with the wing flexure beam fastened to 
a rigid support, because the wing-fuselage junction could take on deformations. 
For comparison, the frequency of each wing mounted on a rigid support and on 
the various flexible fuselages is listed in table l(b). This table includes 
some geometric details of the wing flexure beams and shows the amount of stiff- 
ness reduction due to the flexible wing-fuselage junction. These data are pre- 
sented in figure 10, where it can be seen that the wing-flapping stiffness 
ratio for the wing mounted on the flexible bodies with stiffness levels A and 
B is about 25 and 50 percent lower, respectively, than for the same wing 
mounted on the rigid body. 


Wind- Tunnel-Test Procedure 

The tests were performed in the 2t-inch-diameter leg of the Langley 
Mach 15 hypersonic aeroelasticity tunnel, which uses helium as a test medium. 
This tunnel has a contoured nozzle designed to generate a uniform flow at a 
Mach number of about 15. A description of this facility and its operating 
characteristics can be found in reference 8. 


Models were sting mounted in the tunnel at zero angle of attack. The tun- 
nel with its vacuum reservoir was then evacuated to a pressure of l/2 inch of 
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mercury absolute. A control valve upstream of the test section was opened, and 
flow was established at a dynamic pressure of about 100 lb/sq ft. The dynamic 
pressure was held constant for about 1.5 seconds so that the model transient 
motion could die out. Then with the Mach number remaining constant, the dynamic 
pressure was increased until either flutter was encountered or the maximum tun- 
nel operating conditions were reached. The average total running time for 
these tests was about 5 seconds. 

Throughout the tests, stagnation temperature and pressure were recorded 
on an oscillograph together with signals from resistance-type strain gages 
mounted on the model support system, on the fuselage longitudinal inserts, and 
on the wing flexure beams so that the tunnel conditions could be correlated 
with the model behavior. Flow conditions in the test section were obtained 
from the stagnation pressure and temperature by assuming isentropic nozzle flow 
at the appropriate Mach number. The start of flutter and the flutter frequency 
were determined from the strain-gage response. Fouling switches were mounted 
on the model at appropriate locations to indicate when the amplitude reached a 
value high enough to cause the model to hit the sting. High-speed motion pic- 
tures were taken of the behavior of most of the models during the tests. 

As part of the experimental program concerning this configuration, the 
steady aerodynamic forces and moments were obtained by using an internal strain- 
gage balance mounted in a rigid-body model. A photograph of the model, with 
the upper half of the cylindrical portion of the fuselage removed, and of the 
strain-gage balance can be seen in figure 11. 


METHOD AND APPLICATION OF FLUTTER ANALYSIS 


A Rayleigh-Ritz, coupled-mode flutter analysis, such as that described on 
page 555 of reference was applied to four of the models tested in this study 
and having a wide variation in wing-to-fuselage frequency ratio. The measured- 
mode-shape data for modes 1, 2, 4, and 5 are presented in figures 6 to 9- The 
third mode was not considered because it was antisymmetric and therefore not 
involved in a study of symmetric flutter. 

The determinantal flutter equation used for the flutter calculations Is 
given in the appendix. The integrals in the expressions for the generalized 
mass and aerodynamic coefficients were evaluated numerically by using spanwise 
coordinates 1, 2, 3, and 4 in figure 5> with the chordwise integrands fitted 
to eighth-degree polynomials based on nine equally spaced stations including 
stations on the leading and trailing edges. The generalized mass was based on 
the distribution of masses per unit area calculated from the known weights and 
volumes of the materials used in the construction of the model. The contribu- 
tions from the model support system and flexure beams joining the wing to the 
body were treated additionally as concentrated masses. 

The generalized aerodynamic coefficients given in the appendix are based 
on the quasi- steady aerodynamic approaches of piston theory and Newtonian 
theory. Two applications of these theories were used; one involved the use of 
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piston theory for both wing and fuselage; the other differed from the first 
only in the use of Newtonian theory on the fuselage-nose- cone frustum. Wing- 
fuselage aerodynamic interference effects were not considered. Some calcula- 
tions of the quasi- steady generalized aerodynamic section forces based on 
Newtonian theory as applied to the rounded leading edge of the wing indicated 
negligible contributions of the forces acting on this surface in comparison 
with the piston-theory generalized aerodynamic forces acting on the rest of 
the wing section. As a result of the negligible effect indicated by these cal- 
culations, a separate approximation of the forces acting on the rounded leading 
edge of the wing was not included in the flutter calculations presented in the 
paper. Steady Newtonian aerodynamic forces acting on the spherically blunted 
nose cap of the fuselage were also found to be negligible in comparison with 
those acting on the rest of the nose- cone section. Consequently, a separate 
approximation of the forces acting on the nose cap was also omitted in the 
flutter calculations. 


RESULTS AND DISCUSSION 


Experimental Results 


The results of the wind-tunnel flutter tests are presented in table II, 
which includes a list of the flow conditions at flutter as well as the flutter 
frequency ratio and the flutter-speed-index parameter V for 

each test run. Two of the models tested did not encounter flutter; the data 
given for the models are for the maximum tunnel conditions reached during the 
run. Experimental results from table II are presented in figures 12 and 13 as 
variations of flutter- speed- index parameter and of frequency ratio co/ofy with 


wing-to-fuselage frequency ratio 




The data represented by circles in figure 12 show a flutter speed trend 
for the system in which the frequency ratio has been varied from 0.677 to 1.508, 
while the fuselage stiffness remained constant. This trend indicates that 
increasing the frequency ratio from O.677 to about 1.0 is destabilizing, whereas 
increasing it further, at least to 1.508, is stabilizing. The same general 
trend was obtained when the fuselage longitudinal bending stiffness was varied 
as is shown by the data indicated by the square symbols in figure 12. A rigid- 
fuselage model was tested in order to extend the frequency ratio to zero, and 
the data indicated by the diamond symbols in figure 12 show that its flutter- 
speed-index parameter is almost 60 percent higher than that of the flexible 
model with the lowest frequency ratio (0.677). The flutter- speed- index param- 
eter of the flexible models is always lower than that of the rigid model over 
the range of frequency ratios in the investigation; thus, the flexible model 
never becomes as stable as the rigid- fuselage model. These results show that 
fuselage flexibility is an important flutter parameter and definitely indicate 
the need for simulating fuselage flexibility in flutter studies of this type 
of configuration. 
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The data represented in figure 12 by ticked symbols are for models that 
were retested after the support -system stiffness had been reduced. The sym- 
bols with flags represent data obtained at stiffness level b ( stiffness 
30 percent of original value) and the symbols with tails, data obtained at 
stiffness level c (stiffness 15 percent of original value). The consistency 
of the flutter-speed results, regardless of the support-system stiffness used, 
indicates that the original support-system-stiffness level was low enough to 
simulate the free-flight condition; that is, the support- system rigid-body fre- 
quencies were far enough removed from the frequencies of the elastic modes of 
the model so that they would not influence, or contribute to, the instability 
of the model. 

In figure 13 the flutter frequency as well as the system natural frequen- 
cies, all normalized to the wing flapping frequency, are presented for all 
models tested. The flutter frequency ratio is relatively constant over the 
range of parameters of this investigation. The flutter motion is a coupled 
wing-to-fuselage symmetric mode that changes from being predominantly the 
fourth to predominantly the fifth natural mode for wing-to-fuselage frequency 
ratios greater than about 0.95* 

As part of the experimental program for this configuration the normal 
force, axial force, rolling moment, and pitching moment were measured on a 
rigid- fuselage model at angles of attack of 4° and 10°. These results are pre- 
sented in coefficient form in table III which also includes the aerodynamic- 
center location. Also included in the table are theoretical values for the 
coefficients of normal force and pitching moment about the nose of the model, 
the location of the aerodynamic center, and other pertinent quantities for the 
same configuration. These theoretical forces were determined by using the fol- 
lowing combination of piston theory and Newtonian theory. Newtonian theory was 
applied to the rounded leading edge of the wing, to the spherically blunted nose 
of the fuselage, and to the fuselage cone frustum. Piston theory was applied to 
the planar surfaces of the wing and to the cylindrical portion of the fuselage. 
Shown in table III are the differences between the experimental and theoretical 
forces and moments expressed as percents of the experimental values. These dif- 
ferences are believed to be due to the wing and fuselage aerodynamic forces 
being derived independently without consideration of any wing-fuselage inter- 
ference effects. 


Theoretical Flutter Results and Comparison with Experiment 

The flutter speeds and flutter frequencies calculated by means of the 
coupled-mode flutter equation given in the appendix are plotted in dimensionless 
form along with the experimental flutter speeds and flutter frequencies in fig- 
ures 12 and 13- The flutter calculations were made for the four models whose 
mode shapes are presented in figures 6 to 9* The natural frequencies given in 
these figures for models A-2-a, B-4-a, C-l-a, and A-2-c correspond to those 
given in table I (a) for runs 3 4, 32, 33 > and 37> respectively. Additional flut- 
ter calculations were made for runs 27 , 28, and 31. The mode shapes of 
model A-2-a in figure 6 were assumed to correspond to the natural frequencies 
listed in table l(a) for runs 27 and 28, and the mode shapes of model B-4-a in 
figure 7 were assumed to correspond to the natural frequencies for run 31 * 
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The theoretical results in figures 12 and 13 are based on piston theory 
for the wings and cylindrical fuselage section and on Newtonian theory for 
the fuselage nose-cone frustum. These results were changed by an insignifi- 
cant amount when piston theory was used for all lifting surfaces on the 
configuration. 

From figure 12 it is evident that the theory predicts flutter speed trends 
that are somewhat similar to the experimental trends but are unconservative - 
that is, above experimental trends by about 3 percent to 80 percent. Figure 13 
shows the theoretical flutter frequencies to be in better agreement with experi- 
mental flutter frequencies at wing-to- fuselage frequency ratios near 0.8 than 
for wing-to-fuselage frequency ratios of 0 and above 0.8. It should be noted 
that the values of at which the agreement between theory and experiment 

is poorest are those for which the fifth-mode frequencies, which were the high- 
est measured in this study, approach the experimental flutter frequencies. As 
the spread between these two upper frequencies increases, the theoretical and 
experimental flutter frequencies tend to approach one another. These fre- 
quency comparisons indicate the need for additional higher order modes having 
greater flexibility than that indicated in figures 6 to 9> particularly in the 
chordwise direction on the wing. Unfortunately, the models were destroyed 
during flutter before the flutter calculations were made and before the need 
for more modes was recognized. 

Without the knowledge of additional modes it is impossible to tell whether 
the theoretical flutter-speed trends in figure 12 correspond to converged solu- 
tions, although all the natural modes available were used and found to be 
needed to obtain these trends. For example, in some calculations not presented 
herein, omission of the fourth mode resulted in a 69-percent increase in flut- 
ter speed and a 76-percent decrease in flutter frequency for model A-2-a at 
city 

— — = 0.800 (run 3k ) . Thus, with the limited number of natural modes avail- 
ing. 

able, the adequacies of the quasi-steady aerodynamic approximations involved in 
the flutter mode cannot be truly assessed for the configuration. It can simply 
be noted that such aerodynamic approximations do not account for the actual 
flow field between the strong shock wave and the body surfaces, a fact which 
may well be important. 


CONCLUDING REMARKS 


A wind-tunnel investigation was conducted in helium flow at a Mach number 
of I5A to provide hypersonic flutter data on a low-aspect-ratio, flexible, 
wing-fuselage configuration, and the results are compared with the results of 
flutter calculations. The models were sting mounted in the 24-- inch-diameter 
leg of the Langley Mach 15 hypersonic aeroelasticity tunnel on a support system 
which simulated the free-flight condition by allowing symmetric freedoms of 
vertical, translation and pitch. The fuselage was cylindrical with a spheri- 
cally blunted nose-cone section, and the wing was a 73° clipped delta with a 
modified trailing edge and a slab airfoil section. Model geometry and mass 
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distribution were held constant while the ratio of wing stiffness to fuselage 
stiffness was varied. 

The experimental results indicate a destabilizing flutter-speed trend 
when the wing-to-fuselage frequency ratio is increased from 0 to about 1.0. 

The trend becomes stabilizing when the frequency ratio is increased to values 
greater than 1.0, but the flutter- speed-index parameter of the flexible models 
is always lower than that of the rigid model over the range of wing-to-fuselage 
frequency ratios in this investigation. The results further show that the flut- 
ter frequency ratio is relatively constant over this same range, and that the 
flutter motion is a coupled wing-fuselage symmetric mode. The fact that the 
flutter speed did not appear to be influenced by wide variations in the support 
system stiffness indicates that the original support-system-stiffness level was 
low enough to simulate the symmetric free-flight conditions of vertical trans- 
lation and pitch. The conclusion reached from the experimental flutter results 
is that fuselage flexibility is an important flutter parameter and that flutter 
speeds of a rigid- fuselage model can be unconservative, or higher, than those 
of a flexible-fuselage model. 

Flutter calculations based on the first four measured coupled modes and on 
quasi-steady aerodynamic forces determined by using a combination of piston 
theory and Newtonian theory gave flutter- speed trends similar to experimental 
trends but anywhere from about 3 percent to 80 percent higher. The possibility 
of better correlation between theory and experiment suffers from lack of higher 
order, flexible, wing-fuselage modes with which to approximate the analytical 
flutter mode more accurately and thereby assess the aerodynamic approximation 
more fully. The highly simplified representation of the unsteady aerodynamic 
forces acting on wing and fuselage, with wing-fuselage interference effects 
neglected, is another possible cause of the discrepancy between flutter theory 
and experiment. The use of Newtonian theory for the fuselage nose- cone frustum 
and piston theory for the wing and cylindrical fuselage section resulted in 
essentially the same flutter solution as that obtained by using piston theory 
for all lifting surfaces on the configuration. 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., June 8, 1965. 
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APPENDIX 


COUPLED-MODE FLUTTER EQUATION FOR WING-FUSELAGE CONFIGURATION 
SUBJECTED TO QUASI - STEAD! AERODYNAMIC LOADING 


The theoretical flutter- speed and flutter-frequency trends shown in fig- 
ures 12 and 13 were calculated by solution of the following determinantal 
coupled-mode flutter equation (based on p. 555 of ref. 9 ): 


1 - 
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1 




= 0 (Al) 


1 - 




ii 


CD-. 


where 2 = (l + ig) in terms of reference frequency 0 ) r and structural 


damping coefficient g, and where = 


npb 0 i 2 Z x 

Mi 


in terms of gas density p. 


wing root semi chord b Q , reference length l x , and generalized mass Mj[ which 
is given by 

i = \ mp(x) h i 2 (x) dx + JJ m^(x,y) h i 2 (x,y) dx dy (i = 1, 2, . . .) 

s (A2) 


M. 


The first integral in this equation represents the fuselage contribution 


and the second integral 


= surface integral^, 


the wing contribution to 


the ith generalized mass. The mode shape is denoted by hp; the fuselage mass 
per unit length, by mp; the wing mass per unit area, by iryj and the fuselage 
length, by Ip. 
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APPENDIX 


Aj_j in equation (Al) denotes the general term for the aerodynamic coeffi- 
cients which are also made up of two integrals consistent with the linear super- 
position of forces on the wing and fuselage and the neglect of aerodynamic 
interference effects between the two. With the application of piston theory 
for both wing and fuselage (based on refs. 10 and 11), the general expression 
for A-j^j is 



-§- II h. + i 1* khJ d£ dy (i, j = 1, 2, . . .) (A3) 


dx\Bx 


— _ 

where k = -y-> x = ly = an< ^ X = S X an ^ w ^ ere r is the fuselage 

radius with r = l x r, s is the wing semispan measured from the model center 

line, y is the ratio of specific heats, and t is the local airfoil thick- 
ness. The quantity F (based on ref. 10) is 


= i + 2L±_i M 


(it should be noted that the presence of all positive signs in equation (Ah) in 
contrast with the alternating signs in equation (3*6) in reference 10 is due to 
the difference in sign convention between the present analysis and that of ref- 
erence 10; in the present analysis the coordinate x is considered positive 
rearward in the same direction as the flow.) 

For the present configuration (cone- cylinder fuselage and slab-airfoil- 
section wing) equation (Aj) becomes 




APPENDIX 


where ? = ^x - Zjjjtan 5 + r c in terms of the cone semivertex angle 6, the 
dimensionless fuselage cylinder radius r c , and the dimensionless location of 
the juncture between the fuselage nose cone and cylinder Z N . The reference 
length Z x is chosen as 2b 0 (see fig. l), and all values of x and r are 

referred to this length. Because of the slab airfoil section ^- = 0 and, 

dx 7 

therefore, the last integral in equation (A3) disappears. In the equation for 
F (eq. (A4)) dr/dx is simply tan 8. 

With the introduction of Newtonian theory for the fuselage (ref. 12), the 
first integral in equation (A3) is replaced by 



where 0 is the angle between the local tangent to the surface of the fuselage 
and the longitudinal axis of the fuselage. 

The introduction of Newtonian, theory to the flutter analysis of the pres- 
ent wing-fuselage configuration affects only the first two integrals in equa- 
tion (A5) . With F/M omitted, these integrals are replaced by 



(AT) 

Separate integrations over the blunt leading edges of the fuselage and wing 
were omitted inasmuch as the contributions from these integrals were considered 
insignificant in comparison with contributions from the other integrals. 
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TABLE I.- MODEL PROPERTIES 



Model 

Run 

A-l-a 

29 

A- 2 -a 

34 

A- 3 -a 

27 

A- 4 -a 

28 

A- 5 -a 

30 

B- 4 -a 

32 

B- 5 -a 

31 

A- 2 -b 

36 

A- 2 -c 

37 

C-l-a 

33 

C-l-b 

35 


0.0294 

. 0291 + 

.0294 

.029L 

.0294 

.0294 

.0294 

.0294 

.0294 

.0263 

.0263 






f 3 

f 4 

13.0 

24.0 

12.0 

25.0 

12.5 

27.5 

13.0 

26.8 

13.0 

26.0 

12.5 

20.0 

13.3 

20.0 

9-5 

24.7 

8.0 

24.0 

13.8 

— >00 

12.8 

— >00 




% 

29.3 

102.4 

150.8 

29.0 

125.7 

157.1 

32.0 

150.8 

172.8 

38.5 

207.3 

246.3 

168.4 

163.4 

31.0 

166.5 

125.7 

37.0 

188.5 

125.7 

30.0 

122.5 

155.2 

27.0 

119.4 

150.8 

24.0 

121.9 

00 

23.0 

119.4 

— >00 



0.679 

.800 

.873 

1.231 

1.507 

1.325 

1.500 

.789 

•792 

0 

0 


(b) Wing-f lapping-stiff ness details 



Flexure-beam 

thickness, 

in. 

! 

F 

Model 

Front 

Rear 

f s 

A-l-a 

C-l-a 

C-l-b 

0.012 

.012 

.012 

0.016 

.016 

.016 

20.90 

20.90 

20.90 

A- 2 - a 
A- 2 -b 
A- 2 -c 

o.oi4 

.014 

.014 

0.018 

.018 

.018 

25.00 

25.00 

25.00 

A- 3 -a 

0.016 

0.020 

31.25 

A- 4 -a 

B- 4 -a 

0.020 

.020 

0.025 

.025 

42.00 

42.00 

A- 5 -a 

B- 5 -a 

0.025 

.025 

0.032 

.032 

51.70 

51.70 


Frequencies, cps 


E W,C 1 W,A I W,B 


Stiffness ratios 





16.25 


20.00 

-I FT A 

±.y.jv 

19.00 


24.00 


33.00 


39.20 



0-337 






































































TABLE III.- STEAD! AERODYNAMIC COEFFICIENTS 



Theoretical 
values for 
a = 4° 

(O.O 698 I rad) 

Experimental 
values for - 

Difference between 
theoretical and 
experimental values 
for a = 4°, 
percent of 
experimental value 

a = 4° 

P 

1 ) 

H 

O 

O 

C L 

a o. 00798 




c L a > per rad 

0.1143 




q., lb/sq. in 

4.02 




Z M 2 , in . 2 

441 




l?> in 

18 




C A 


0.0030 

0.0034 


c z 


0 

0.0010 


Cn 

0.0080 

0.0086 

0.0248 

-7.0 


0.00427 

0.00514 

0.01528 

-17.0 

x ac 

0.534 

0.597 

0.617 

-10.6 
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(a) Cone and integral longitudinal inserts with unmounted coil spring. L-64-3338 
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(c) Side view of model with wings mounted and rubber covering partially installed on fuselage. L-64-3335 




2k 


Figure 2.- Continued. 
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(e) View of wing showing root chord and trailing edge. L-64-3336. 1 
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(f) Rigid-fuselage model with upper portion of cylindrical fuselage rotated out of place. L-64-3337 

Figure 2.- Concluded. 







Weight of fuselage per unit length, Ib/in. 



Wing — 

Concentrated mass due to flexure-beam weight, lb, at- 

flapping- 

X = 8.5 in. X = 17.5 in. 

stiffness 

0 A O D O A 

level 


1 

0.00025 0.00051 0.00087 0.00087 0.00087 0.00087 

3 

.00034 .00067 .00110 .00108 .00108 .00108 

4 

.00042 .00064 .00135 .00135 .00135 .00135 

5 

00052 .00102 .00173 .00173 .00173 .00173 


Spanwise 
station, in. 
-0-1.5 
-a- 2.5 
-0-3.5 
— A-4.5 



6 10 14 18 2 

Chordwise station, in. 


(c) One wing. 


Figure 4.- Concluded. 




1234 565 4321 


Spanwise coordinate, y , in. 


Figure 5 - Coordinate system and model center-of-gravity location. 
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coordinate 
o | 

□ 2 
O 3 
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Q 6 



4 8 12 16 


Chordwise coordinate 


(e) Fifth mode (wing-fuselage combination); f$ 


Figure 6.- Concluded. 



Displacement normalized to maximum value in a given mode 


Spanw i se 
coordinate 

Q 6 



0 4 8 12 16 


(a) First mode (model translation); ^ = 7.10 cps. 



0 4 8 12 16 


Chordwise coordinate 


(b; Second mode (mode! pitch); f 2 = 10.7 cps. 
Figure 7.- Mode shapes for model B-4-a. 



(c) Fourth mode (fuselage bending); = 20.0 





4 8 12 16 

Chordwise coordinate 


(d) Uncoupled wing flapping mode (fuselage restrained); l 


Figure 7.- Continued. 




Spanwise 
coordinate 
o I 
a 4 
o 6 



4 8 12 16 


Chordwise coordinate 


(e) Fifth mode (wing-fuselage combination); f 5 = 31.00 cps. 
Figure 7.- Concluded. 





Spanwise 
coordinate 
o l 
a 2 
<0 3 
a 4 



4 8 12 16 


(c) Uncoupled wing flapping mode (fuselage restrained); fyy = 19.^ 



cr-' 0 ' o 

4 8 12 16 
Chordwise coordinate 


id) Fifth mode (wing fuselage combination); f 5 = ?4.0f)cps. 


Figure 8.- Concluded. 










Fuselage 

stiffness 

level 

o C 



2 3 4 5 

Wing-flapping-stiffness level 


43 


Figure 10.- Reduction in wing-flapping-stiffness ratio due to flexible wing-fuselage junction for various models. 



Figure 11.- Aerodynamic-force model and strain-gage balance. 


L-63-3789. 1 




Fuseloge 
stiffness level 
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O C 
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Typicol Meaning 
symbol 

• No flutter 

d Support system 

stiffness b 

q Support system 

stiffness c 


i i i 1 1 

.4 .8 1.2 1.6 2.0 

Wing- to- fuselage frequency ratio, ^ W 


Figure 12.- Variation of flutter speed with wing-fuselage frequency ratio. 



